Magnetic field as a probe of gap energy scales in YBCO 
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Among high Tc materials, the YBCO (YBa2Cu307_a;) compounds are special since they have 
superconducting chains as well as planes. We show that a discontinuity in the density of states as a 
function of magnetic field may appear at a new energy scale, characteristic of the chain and distinct 
from that set by the d wave gap. This is observable in experimental studies of the thermodynamical 
properties of these systems, such as the specific heat. 



It is widely accepted that a d-wave pairing gap de- 
scribes the superconductivity in high Tc compounds [Q. 
The YBCO compounds are special however, due to the 
fact that the one dimensional chains change the Fermi 
surface, giving a quasi-linear branch, and thereby intro- 
ducing a new energy scale for the variation of the gap 
on the Fermi surface. We argue that this can be directly 
measured using a magnetic field. 

Scanning tunnelling microscopy (STM) measurements 
[ [H,D of the local density of states of YBCO agree with 
earlier tunnelling measurements [Q and show a very clear 
structure around 20 meV corresponding to the full gap 
energy. They also indicate structure at around 5.5 meV 
which is absent in BSCCO compounds j^]. It is natural 
to assume that this is because YBCO has a different elec- 
tronic structure, possessing orthorhombic chains as well 
as tetragonal planes. The chains introduce a new energy 
scale, the value of the gap function at the endpoint of 
the chain Fermi surface on the boundary of the Brillouin 
zone. Since experiments are only sensitive to the value 
of the gap function on the Fermi surface (and not to its 
global behaviour throughout the Brillouin zone) this en- 
ergy scale will be apparent in experiments. Another in- 
terpretation of the additional structure is localised quasi- 
particle states in the vortex core ||]. Even if present, 
these should have a negligible effect on the bulk density 
of states, a point also made in Ref. 

We model this by assigning the chains an elliptical 
Fermi surface but with a highly modified gap function 



Ao(cos 2(f> + s). 
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(j) is the polar coordinate in k space and the dimension- 
less parameter s (with \s\ < 1) indicates the effective 
amount of s wave in the order parameter. We stress that 
(|l|) is not the gap function defined throughout the Bril- 
louin zone but is rather its value evaluated on the Fermi 
surface. Along the chain Fermi surface the gap function 
will have a dependence like (|l|) even if it is globally pure 
d-wave. This is an effective theory to model the contri- 
bution of the chains to the density of states and reflects 
a combination of gap symmetry and Fermi surface ge- 
ometry. (Independent information on the existence of a 
subdominant s-wave component is provided by Joseph- 



son tunnelling data 0,^.) 

Standard theory for the density of states gives 
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where is the dispersion relation for the quasi-particle 
excitation energy (for the moment we take it = k^/2m 
with k the momentum and m the effective mass) and N is 
the normal, nonsuperconducting density of states. Using 
the gap function (^), this integral is easily done and is 
plotted in Fig. ^ using s — 0.57 (this choice is explained 
below.) There is a van Hove singularity at a; = (l + |s|)Ao 
as we expect ||l^], however there is another at (1— |s|)Ao, 
in qualitative agreement with 




FIG. 1. Typical result for N{uj) for a pure superconductor 
using s = 0.57 to introduce a second low energy scale. Note 
the van Hove singularities at oj/Aq = (1 ± |s|). 

A problem with the measured N{u!) is that there is a 
large zero bias anomaly in contradiction with theory and 
presumably arising from surface effects. It is therefore 
useful and possibly better to measure a bulk property. 
We propose the low temperature limit of the specific heat, 
which comes from the density of quasiparticle states at 
zero frequency. We focus on the magnetic field regime 
Hci <^ H Hc2 where the density of states is given by 
quasi-particle excitations in the presence of the induced 
vortices. For each magnetic field, there is a magnetic 
energy, which scales as ^/H [|11|,^ . We expect structure 
at critical fields for which the magnetic energy sweeps 
through the singularities at Ao(l± s) in Fig. |^. The new, 
lower energy scale is well within the range of available 
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fields and is observable in specific heat experiments which 
probe the vortex modified density of states. 

In the presence of the vortex, the Cooper pairs have 
a superfluid velocity Vg which depends inversely on the 
distance to the vortex core. The quasi-particle energies 
of momentum k are Doppler shifted by a; — > + Vg • k 
I p^ , which can be thought of as introducing a spatially 
dependent shift in chemical potential. We are interested 
in the zero frequency case, so we can use Eq. (||) but 
with u; replaced by y = Vg • k. In addition to the trace 
integral over k, we also do a spatial average within the 
vortex unit cell, with coordinates r: 
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To find Vg, we use the free energy density per unit 
length 



F = / d 
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We have introduced the parameters A| = /i^A^ with 
= Mc^/ine^Us, fii = nii/M and M = ^nixmy (A 
is a mean London penetration depth and Ug the super- 
fluid density.) Following standard free energy minimisa- 
tion fl^ with a flux line source at the origin, we express 
the magnetic field as a modified Bessel function. In the 
magnetic field regime considered here, the typical spac- 
ing between vortices is much smaller than A and we use 
the small argument approximation: 
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pointing in the z direction and $0 — ixhcje is the flux 
quantum. Applying Ampere's law, the current is 
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We find Vg by dividing the current by eng. 

The quasi-particle excitations are given by the disper- 
sion relation 
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As discussed above, we take the gap function (|l|) to be 

(8) 



Ak = Ao/(k) = Ao I S— S + s 



The analysis is simplified by a change of variables to 
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In terms of these new coordinates we have 
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where {r',(3') are the spatial polar coordinates {(3' is 
the vortex winding angle) and {k',(f)') are the momen- 
tum polar coordinates in the new coordinates. Note that 
Vg • k = |vg||k| sin(/3' — cj)'). Henceforth we drop the 
primes. Since we are integrating over /3, we are also free 
to shift its origin and thereby replace f3 — cj) hy [3. To 
conform to usual notation we say that the chains run 
along the h direction which we define as parallel to the 
y axis. Since the carrier mass in the h direction is less 
than in the a direction (due to the conductivity supplied 
by the chains) we conclude that niy/mx < 1. In prin- 
ciple s can have either sign (we ignore the time-reversal 
symmetry breaking possibility for which s is complex.) 
For purposes of exposition we will use the realistic value 
'niy/mx = 0.5 and take s = 0.57 so as to get a second gap 
scale around 5.5 meV as seen in experiment. The oscil- 
latory part of the gap function in ( px| ) has been modified 
due to the coordinate rescaling, compared to the origi- 
nal Aq cos 2(j). Nevertheless, the amplitude has not been 
modified and without the s term both functions would 
vary between plus and minus Aq. Either form would 
yield an initial linear dependence on 12 1 followed 

by a saturation at a scale set by Aq. Rather, it is the s 
term which changes this qualitative picture. It makes the 
gap function have a maximum value of Ao(l + |s|) and a 
minimum value of Aq (1 — |s|). The first of these gives the 
saturation scale mentioned above while the second gives 
structure corresponding to the first peak in Fig. 

We begin by nondimcnsionalising the spatial integral 
in (H). In order for the average magnetic field in the sam- 
ple to equal the applied field, we require the inter-vortex 
spacing to be [|l2| R — ^J<^q/iiH /a where a is a geo- 
metrical constant of order unity which accounts for the 
mismatch between the circular vortices and the hexago- 
nal lattice they fill out. We define the magnetic energy 
as 
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where vp is the Fermi velocity and the dimensionless 
magnetic energy as = Eh/A.q. We then express (||) 
as 
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where N is the normal density of states and p — r/R. 
We used the delta function to perform the k integral and 
it is understood that the integration domain is limited to 
the range where the integrand is real. 

We can do the p and (3 integrals in ( ]T^ ) leading to the 
final expression for the density of states, 
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The parameter v is experimentally tunable while s (which 
enters f{4>) through (10)) is a fixed, intrinsic material 
parameter (although it may depend somewhat on doping 
levels) and represents a combination of gap and Fermi 
surface symmetries. For the case |s| < 1 and v <^ 1, it 
is simpler to go back to ( p^ ) and expand the integrand 
around the gap nodes. One can perform the integral 
approximately in the limit of small v and the resulting /3 
and p integrals are trivial so that 



(14) 



The index n refers to the nodes of the gap function and 
/4 is the derivative of the function / evaluated at the 
nodes. This is still linear, as one has for a pure d-wave 
gap, although the presence of the s term does change the 
positions of the gap nodes and hence affects the slope. 
For \s\ > 1, there are no nodes and the dependence is 
quadratic for small v. Since this is unphysical we do not 
present the explicit form. 
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FIG. 2. The four possible behaviours of the integrand of 
Eq. ( |l3| ) as a function of (p. The values of {u, s) are (0.4, 0.57), 
(0.5,0.57), (2.0,0.57) and (0.5,1.6) for cases a to d respec- 
tively. In all cases the mass anisotropy was my/m^ = 0.5. 
Top: regions in the v — s plane corresponding to the 4 types 
of integrands (in regions b' and d' the integrand is translated 
by 7r/2 from its form in regions b and d.) 

There are four distinct behaviours of the integrand of 
( [Tsl ) , as shown in Fig. || and the v — s plane is accord- 
ingly divided into six regions. For fixed s, we vary v and 
thereby cross from one region to another. Associated 
with such crossings there is a corresponding nonanalytic- 
ity in the density of states, which is experimentally mea- 
surable. This is shown in Fig. ^ where we plot the density 
of states as a function of v using the integral form (p^). 
Clearly there is structure ai v = 0.43 and 1.57, in accord 
with our general considerations. The kink a,i v = 1.57 
corresponds to the integrand saturating at unity — as 
we go from an integrand as in Fig. ^3 to the type as in 
Fig. However, this probably lies beyond any realis- 
able magnetic field strength. The kink at v — 0.43 which 
corresponds to going from an integrand as in Fig. ^ to 
the type in Fig. appears by eye to be a discontinuity 
in slope. This is not really the case, as we now discuss. 

The integrand in Fig. ^ja saturates around (j) — t: /2 and 
(j) = 37r/2 when v approaches 1 — |s| and the integral does 
not increase as quickly after this value as it does before. 
We can evaluate the missing contribution, 5N , from the 
regions around these local minima. We take u to be just 
below this transition, i.e. 



V = 1 



(15) 



with e <C 1. To leading order in e, the contribution of 
this region to (|l^) is 
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There is no such contribution if v is just above 1 — |s|, so 
there is a discontinuity in the density of states. While the 
precise prefactor and position of the discontinuity may 
depend on factors which we have not included, the three- 
halves power is generic, depending only on the topological 
property of having local minima of the integrand disap- 
pear. For example, in the event that s < 0, we first lose 
the minimum around (j) = Q instead of (j) — tt/2 but oth- 
erwise the behaviour is the same and, except for a change 
in prefactor (typically it is smaller), we still expect a dis- 
continuity of this power. Similarly, it is the same order 
of discontinuity at = 1 -I- |s| although, as stated, this is 
probably beyond the accessible magnetic field range. 

To compare with experiment, we now determine the 
parameters of our model. The upper scale of 20 meV 
equals (1 + s)Ao and the lower scale of 5.5 meV equals 
(1 — s)Ao from which we infer Aq ~ 12.7meV and s « 
0.57. We also use the result, consistent with both /uSR 
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and specific heat measurements that Eh ~ 2.6 vff when 
H is measured in Tesla and Eh in meV JTH . 

The dimensionless parameter ly then equals 0.20%///. 
The dependence of the specific heat on H has been mea- 
sured in p5|-[l^, with having the cleanest samples. 
The data is shown in FigTH in dimensionless units (taking 
the specific heat at low temperatures to be proportional 
to temperature and to the density of states.) By our 
previous arguments, energies of 5.5 meV correspond to 
magnetic fields of about 4.5T or of about 0.42, which 
is well within the experimental ranges considered. The 
origin of the finite value of the density of states at zero 
field is not well understood but is sometimes attributed 
to oxygen vacancies on the chain [p^ or to resonant im- 
purity scattering dj. Either way, to compare with our 
results one should subtract off this constant. There will 
also be a linearly increasing component from the den- 
sity of states on the planes; this should be of comparable 
magnitude but without the lower singularity. 




0.4 



5^ 0.2 ■ 



0.0 



0.0 



0.4 0.6 
V 



1.0 



FIG. 3. Top: the density of states plotted against the 
rescaled magnetic field v for the case s = 0.57 and 
my/rrix = 0.5. Bottom: experimental data converted to di- 
mensionless units as described in the text. Circles, diamonds 
and squares are from [15,16,18] respectively. In each case we 
have estimated a best-fit line below and above f = 0.5 to 
guide the eye. 

We plot this experimental data, fitting it with straight 
lines of different slopes at the two extremes of the data 
sets. There does appear to be a change in slope around 
ly — 0.5, in qualitative agreement with theory. While the 
theory does not predict a discontinuous slope, it does pre- 
dict a nonanalyticity which resembles a change of slope. 
Since the data points are so sparse, we have not tried 
a detailed fitting to the predicted functional form. The 
magnitude of this nonanalyticity is given by the prefac- 
tor of the e'^/^ term in ( p^ . (We note that the aforemen- 
tioned background slope from the planes is present in the 
experimental data but not in the theoretical curve, so the 
relative magnitude of the discontinuity is different.) Ex- 
periments using a much finer sampling of magnetic field 
values would be required to verify this prediction. 

In this paper we have considered a magnetic field paral- 
lel to the c-axis. In the case where it is parallel to the a—b 



planes there may be additional interesting anisotropy ef- 
fects [0. Unfortunately, the energy scales seem to be 
such as to put this out of the experimentally accessi- 
ble magnetic field range |2^. An interesting possibil- 
ity is to consider a tilted magnetic field such that the 
energy scales are experimentally accessible while the in- 
plane component is still strong enough to yield observable 
anisotropy effects. A gap function of mixed symmetry 
would have a clear signature on this anisotropy. Another 
interesting extension is to consider the role of the param- 
agnetic response of the electrons to the magnetic field Q 
— an effect which has been completely neglected in the 
present work. While the energy scales seem to be such 
as to make this a reasonable assumption, it could well be 
that for the in-plane magnetic field, the paramagnetism 
is of comparable importance. 
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